Conformational phase diagram for polymers adsorbed at ultrathin nanowires 
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We study the conformational behavior of a polymer adsorbed at an attractive nanostring and 
construct the complete structural phase diagram in dependence of the binding strength and effec- 
tive thickness of the string. For this purpose, Monte Carlo optimization techniques are employed 
to identify lowest-energy structures for a coarse-grained hybrid polymer-wire model. Among the 
representative conformations in the different phases are, for example, compact droplets attached 
to the string and also nanotube-like monolayer films wrapping the string in a very ordered way. 
We here systematically analyze low-energy shapes and structural order parameters to elucidate the 
transitions between the structural phases. 

PACS numbers: 82.35.Gh,05.10.Ln,61.48.De 



The interaction of polymers with substrates is relevant 
in nanoscale applications like molecular nanoelectronic 
circuits and in biological processes such as receptor- 
ligand binding. Thus, scrutinizing basic structural mech- 
anisms of molecular binding at interfaces is crucial for a 
large field of interdisciplinary research and for potential 
applications. 

In recent years, there has been substantial progress in 
understanding general properties of polymer adhesion at 
solid substrates. This includes, for example, the identi- 
fication of generic structural phases and the transitions 
between these [H-Q, as well as specific binding affinities 
of proteins regarding the t ype of the substrate and the 
amino acid sequence (il. Ho 12]. In most of these studies, 
the substrate is considered as being planar. The influ- 
ence of curved substrates on the formation of structural 
phases has been subject of works on droplet and helix 
formation at cylinders [HI, [l4| . 

In this Letter, we systematically study conformational 
phases induced by an attractive nanowire, i.e., a sub- 
strate with one-dimensional topology. A nanowire could 
be, for example, a stretched polymer with the ends at- 
tached to dielectrical beads fixed by optical tweezers. It 
is one of the most striking results of our study that under 
certain parametrizations of the polymer-nanowire inter- 
action, i.e., in the corresponding region of the conforma- 
tional phase diagram, the polymer crystallizes in stable 
cylindrical shapes with monomer alignments which re- 
semble atomic arrangements known from single-walled 
carbon nanotubes. In this conformational phase, the 
cylindrical hull surrounds the thin wire such that the in- 
terior is free of particles. In a hydrodynamic application, 
for example, molecules can still flow through it. Since 
the axis of a polymeric tube is always oriented parallel 
to the direction of the wire, the growth direction of the 
tube can be controlled. This would enable the construc- 
tion of complex tube systems and, therefore, allows for 
applications beyond those known for conventional nano- 
tubes. Another conceivable application of this structural 
coincidence is the systematic stabilization or functional- 



ization of nanotubes by polymer coating 15-17). 

In our study, we investigate a coarse-grained model 
for the polymer and a linelike substrate representing the 
nanowire. For the polymer, we employ a linear bead- stick 
model, i.e., covalent bonds between the N monomers are 
stiff. The chain is not grafted to this string and may 
move freely. The total energy of the system includes 
three contributions, E = + £?bend + ^string- The 
interaction between nonadjacent monomers is governed 
by the standard Lennard- Jones (LJ) potential, 
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with the distance between nonbonded monomers i 
and j. The monomer-monomer interaction parameters 
e m and <r m are set to unity in the following. The weak 
bending energy is a remnant of the protein-like origin of 
the model [l9[ and reads ^bend({cos^}) = i^Y^i^i^ ~ 
cos 6i) with the bending stiffness set to k = 1/4. The 
bending angle Oi is defined by the covalent bonds con- 
nected to the ith monomer. The monomer-string energy 
is obtained by continuously integrating a standard LJ 
potential over the infinitely long string [18]. We find 



N 



£string({r± ; i}) =7Taef^2 




(2) 



where <jf and ef are the monomer-string interaction pa- 
rameters. The distance of the ith monomer perpendicu- 
lar to the string is denoted by r± ] i. For convenience, we 
scale the potential such that its minimum is —1 at r™ 
for ef = 1 and <jf = 1, in which case a ~ 0.528 [18]. The 
effective thickness of the string, erf, is related to the min- 
imum distance r™ of the monomer-string potential via 
r^ in (cr f ) = (693/480) 1/6 a f « 1.06 a f . Alternatively, the 
monomer-string energy can be considered as the limit- 
ing case of the interaction of a monomer with a cylinder 
(cp. Ref. [Hj) of radius R —> 0, keeping the overall LJ 
"charge" fixed 0. 
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FIG. 1: The conformational phase diagram parametrized by 
the monomer-string potential parameters; from left to right, 
the effective string thickness <jf increases and from bottom 
to top, the string attraction strength ef gets larger. Gray 
bands (widths correspond to uncertainty) indicate transition 
lines between compact, crystalline polymer structures with 
the string inclosed (Gi) or excluded (Ge), crescent-shaped 
(C), and barrel- like (B) conformations. The dashed line indi- 
cates a topological crossover that separates mono- and mul- 
tilayer regions. Inset pictures show representative low-energy 
states. Monomers with the same coloring (or shadings) belong 
to the same layer. 



We now systematically analyze the conformational 
phases of a polymer with N = 100 monomers for differ- 
ent values of effective thickness and attraction strength of 
the string, <jf and ef , respectively. By "phase" we denote 
a domain in the parameter space, where the representa- 
tive conformations share qualitatively the same morphol- 
ogy. We have convinced ourselves in exemplified simula- 
tions that the results we discuss in the following for the 
lOOmer are also qualitatively correct for longer chains. 
In Fig. [TJ the conformational phase-diagram is shown 
and representative adsorbed polymer structures are de- 
picted. This phase diagram is a result of extensive analy- 
ses of structural properties for more than 150 low-energy 
conformations with different parametrizations. For the 
identification of lowest-energy conformations, stochastic 
generalized-ensemble Monte Carlo methods [20| and de- 
terministic conjugate-gradient optimization were used. 
Conformational changes were performed by applying 
a variety of update moves, including local crankshaft, 
slithering-snake, global spherical-cap, and translation 
moves [18j . 

In the phase diagram, four major structural phases 
can be identified. For weak attraction (ef < 3), two 
types of crystalline droplets [21[ adhered to the string 
(regions Ge, Gi) can clearly be distinguished. Either 
the string axis is inclosed inside the droplet (Gi) or 
passes by externally (Ge). If the adhesion strength of 
the string ef increases, compact droplets in Gi melt near 
ef « 3 and phase B is entered, where polymer confor- 
mations extend along the string axis. Near ef ~ 4.5, 




FIG. 2: Asymmetry parameter A at small effective string 
thickness (<jf = 0.5). The inset pictures show correspond- 
ing conformations at ef = 2.0, 3.5, and 4.5, illustrating the 
transition from spherical-symmetric to stretched cylindrical- 
symmetric low-energy structures. The conformational tran- 
sition between Gi and B occurs near ef = 3.0. For larger 
values of ef, A starts to significantly deviate from zero and 
conformations become cylindrical. 



a crossover (dashed line in Fig. [T]) from the multilayer 
barrel structures to monolayer conformations with strong 
similarities to single- walled nanotubes occur. According 
to former studies of polymer adsorption at planar sub- 
strates 0, [s[ , this crossover corresponds to a topological 
transition between three-dimensional compact crystalline 
and two-dimensional filmlike structures. For sufficiently 
large values of the effective string thickness erf, the poly- 
mer layers do not completely wrap the tube and stable 
crescent-shaped "clamshell-like" [2[ structures dominate 
in the region denoted by C. 

Let us now have a closer look at the different confor- 
mational transitions. In region Gi, compact crystalline 
conformations with spherical symmetry dominate. In- 
creasing in this regime the attraction strength ef while 
keeping the effective thickness erf fixed, this symmetry 
breaks at ef ~ 3 and the cylindrical phase B is entered. 
This transition can be best characterized by introduc- 
ing an asymmetry parameter based on the gyration ten- 
sor components parallel and perpendicular to the string, 
A = r| yr /r^ yr — 1. This order parameter is shown as a 
function of ef in Fig. [2j exemplified for df = 0.5. In the 
spherical regime Gi, A « 0. As expected, A increases for 
ef > 3 and the structures become asymmetric. At this 
point, it is equally favorable for a monomer to stick to the 
string, or to form contacts to neighboring monomer lay- 
ers. The conformations stretch along the string until they 
form a maximally compact monolayer tube surrounding 
the string for ef > 4.5. 

Increasing, on the other hand, the effective thickness erf 
for values of the attraction strength ef < 3, the transition 
from Gi to Ge is characterized by the different locations 
of the string relative to the droplet: it is inclosed (Gi) 



3 



£f 

0.5 1.0 1.5 2.0 2.5 3.0 




I I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ L 

0.5 0.6 0.7 0.8 

Of 



FIG. 3: The plot shows the opening angle a of low-energy 
conformations for ef = 3/2 in dependence of <jf (squares, lower 
scale) and for erf = 5/6 as a function of ef (circles, upper 
scale). The inset pictures show corresponding conformations 
at ef = 3/2 and <jf = 1/2,4/5, illustrating the separation of 
the droplet from the string. 
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FIG. 4: Distance Ar of the center of mass of the polymer 
from the virtual surface of the cylinder with the radius that 
corresponds to the minimum position of the string potential 
(r™ m ~ 1.06crf) for erf = 7/3. The intersection of the curve 
with the dotted line (Ar = 0) at ef = 2.9, where the center of 
mass equals the radius of this cylinder, defines the transition 
from Ge to C. Pictures show conformations at ef = 2, 3, and 4. 



or excluded (Ge). For small values of ef, the transition 
point can be estimated as a first approximation by as- 
suming a tetrahedral monomer-packing in the crystalline 
droplet. Then, the circumsphere radius of a tetrahedron 
is r Q ~ 0.61 which corresponds to a limiting effective 
string thickness <7f j0 « 0.58. Thus, inserting a string 
with erf < <jf j0 does not break intra-monomer contacts 
within a compact structure. Above this limiting value, 
however, the string would cause an energetically disfa- 
vored replacement of monomers inside the conformation 
and is hence "pushed" out of the droplet. 

Quantitatively, this transition can be identified by 
measuring the opening angle a of a given conforma- 
tion. Projecting the positions of monomers in contact 
with the string onto a plane perpendicular to the string, 
a is defined as the angle between the string and two 
monomers that spans the largest region of the plane with 
no monomers residing in. Thus, roughly, conformations 
with a < 7r correspond to conformations inclosing the 
string (Gi), whereas a > if the string is located out- 
side the droplet (Ge). Figure [3] shows how a changes 
when crossing the transition line Gi^Ge horizontally or 
vertically. Fixing the attraction strength at ef = 3/2, a 
increases rapidly from erf ~ 0.65 (squares, lower scale), 
which is close to the estimate <jf >0 given above. For larger 
values of the effective thickness the string is shifted out- 
wards to retain optimal monomer packing. The inset pic- 
tures show lowest-energy conformations at ef = 3/2 and 
cr f = 1/2 (representative for phase Gi) and erf = 4/5 (Ge), 
respectively. Increasing, on the other hand, the string at- 
traction strength ef while erf = 5/6 is fixed (circles, upper 
scale), the inclusion of the string, accompanied by a rapid 
decrease of a, occurs at ef ~ 1.75. 

Starting in Ge and increasing erf and ef above cer- 
tain threshold values results in the transition towards 



adsorbed curved conformations (C) in the sense that the 
polymer begins to wrap the string. Different monomer 
layers form. We quantitatively define this transition 
to occur at the point, where the distance of the cen- 
ter of mass of the polymer from the string, r^° m = 
Y^iLi e( l ua l s the monomer-string potential 

minimum distance r™ m (~ 1.06<7f), i.e., at Ar = r c f m — 
r mm _ q Qualitatively, the center of mass intrudes into 
the virtual cylinder with radius r™ m , defined by the in- 
ner layer of monomers. In Fig. 2J Ar is plotted as a 
function of ef at erf = 7/3. The transition point Ar = 
is marked by the dotted line which is here intersected at 
ef = 2.9, in correspondence to the Ge<H>C transition line 
in the phase diagram in Fig. [TJ The shown inset pic- 
tures represent conformations with Ar = 0.6,-0.1,-0.8 
at e f = 2,3,4. 

Finally, increasing ef further, region B is entered, 
i.e., ground-state polymer conformations wrap the string 
completely. If the attraction between a monomer and 
the string becomes stronger than the interaction between 
stacked, neighboring monomer layers, regular monolayer 
films surrounding the string are formed, i.e., single- walled 
tubes with an ordered arrangement of monomers. It is 
noticeable that there is a competition between different 
chiral orientations of the wrapping in dependence of the 
monomer-string interaction length scale erf [18]. This be- 
havior is in a similar manner known from carbon nano- 
tubes [22]. Defining the wrapping vector c= nti\ + ma2, 
with di and <22 being the unit vectors of the structure, 
we even find the limiting "armchair" and "zigzag" struc- 
tures, corresponding to m = n and m = 0. Examples 
of nanotube-like polymer conformations with different 
chiralities are shown for erf = 1.50, 1.57 in Figs. [5ja) 
and EJb), respectively, and for erf = 0.65 in Fig. [5fc) 
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FIG. 5: Highly ordered cylindrical monolayer conformations 
of adsorbed polymers with different wrappings in the barrel 
phase B at ef = 5 for (a) <jf = 1.50, (b) <jf = 1.57, and 
(c) erf = 0.65 (different colors or shadings shall facilitate the 
perception only). Geometric properties of these structures 
resemble chiral alignments of atomic structures known from 
single- walled nanotubes. 

(all at ef = 5). The alignment of monomers in Fig. [5ja) 
is almost parallel to the string, whereas the conforma- 
tion in Fig. Efb) exhibits a noticeable chiral winding. If 
the radius of the monolayer polymer tube does not al- 
low for a perfect monomer alignment, defects occur and 
cause the formation of structural domains with different 
chiralities within the same conformation [18|, as in the 
example shown in Fig. [5fc). 

To summarize, we have constructed the entire confor- 
mational phase diagram of a hybrid system consisting of 
a flexible polymer and an ultrathin attractive nanowire 
in dependence of the energy scales and length scales asso- 
ciated to the polymer-nanowire interaction. We identi- 
fied conformational phases of compact spherical polymer 
droplets inclosing or excluding the string, and a phase of 
compact but curved shapes (crescent-shaped structures). 
For sufficiently large string attraction strengths, we ob- 
serve the formation of cylindrical conformations which in 
the extreme case of monolayer structures possess strong 
similarities to nanotubes. This is particularly interesting 
as it shows that polymers can form tubelike structures 
in a controlled way. Since the polymer tube can adapt 
any orientation of the guiding nanowire, also the forma- 
tion of complex, nonlinear tube systems with bends is 
conceivable. This would enable a wide range of poten- 
tial applications which are hard to construct by atomic 
nanotubes. 
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